Long-range interaction and nonlinear localized modes in photonic crystal waveguides 
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We develop the theory of nonlinear localised modes {intrinsic localised modes or discrete breathers) 
in two-dimensional (2D) photonic crystal waveguides. We consider different geometries of the waveg- 
uides created by an array of nonlinear dielectric rods in an otherwise perfect linear 2D photonic crys- 
tal, and demonstrate that the effective interaction in such waveguides is nonlocal, being described 
by a new type of nonlinear lattice models with long-range coupling and nonlocal nonlinearity. We 
reveal the existence of different types of nonlinear guided modes which are also localised in the 
waveguide direction, and describe their unique properties including bistability. 
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I. INTRODUCTION 



In physics, the idea of localisation is generally associ- 
ated with disorder that breaks translational invariance. 
However, research in recent years has demonstrated that 
localisation can occur in the absence of any disorder and 
solely due to nonlinearity, in the form of ^intrinsic lo- 
calised modes, also called discrete breathersu A rigorous 
proof of the existence of time-periodic, spatially localised 
solutions describing such nonlinear modes has been pre- 
sented for a broadxlass of Hamiltonian coupled-oscillator 
nonlinear lattices,B but approximate analytical solutions 
can also be found in many other cases, demonstrating a 
generality of the concept of nonlinear localised modes. 

Nonlinear localised modes can be easily identified in 
numerical molecular-dynamics simulations in many dif- 
ferent physical models (see, e.g., Ref. |l| for a review), 
but only very recently the first experimental observa- 
tions of spatially localised nonlinear modes have bee: 
reported in mixed-valence transition metal complexes,! 
quasi-one-dimensional antifciftomagnetic chains^ and ar- 
rays of Josephson junctions.EI Importantly, very similar 
types of spatially localised nonlinear modes haver-been 
experimentally observed in macroscopic mechanicalEl and 
guided-wave optical] systems. 

Recent experimental observations of nonlinear lo- 
calised modes, as well as numerous theoretical results, 
indicate that both effects, i.e. nonlinearity-induccd lo- 
calisation and spatially localised modes, can be expected 
in physical systems of very different nature. From the 
viewpoint of possible practical applications, self- localised 
states in optics seem to be the most promising ones; they 
can lead to different types of nonlinear all-optical switch- 
ing devices where light manipulates and controls light 
itself, by varying the input intensity. As a result, the 
study of nonlinear localised modes in photonic structures 
is expected to bring a variety of realistic applications of 
intrinsic localised modes. 

One of the promising fields where the concept of non- 
linear localised modes may find practical applications is 
in the physics of photonic crystals [or photonic band gap 



(PBG) materials] — periodic dielectric structures that 
produce many of the same phenomena for photcas as the 
crystalline atomic potential does for electrons^ Three- 
dimensional (3D) photonic crystals for visible light have 
been successfully fabricated only within the past year or 
two, and presently many research groups are working on 
creating tunable band-gap switches and transistors op- 
crating entirely with light. The most recent idea is to 
employ nonlinear properties of band-gap materials, thus 
creating nonlinear photonic crystals that have 2D or 3D 
periodic nonlinear susceptibilityMLj 

Nonlinear photonic crystals or photonic crystals with 
embedded nonlinear impurities create an ideal environ- 
ment for the generation and observation of nonlinear 
localised photonic modes. In particular, such modes 
can be excited at nonlinear interfaces with quadratic 
nonlinearityjij or along dielectric waveguide structures 
possessing a nonlinear Kerr-type response.EHI In this pa- 
per, we analyse nonlinear localised modes in 2D photonic 
crystal waveguides. We consider the waveguides created 
by an array of dielectric rods in an otherwise perfect 2D 
photonic crystal. It is assumed that the dielectric con- 
stant of the waveguide rods depends on the field intensity 
(due to the Kerr effect), so that waveguides of differ- 
ent geometries can support a variety of nonlinear guided 
modes. We demonstrate here that localisation can oc- 
cur in the propagation direction creating a 2D spatially 
localised mode (see Fig. || below). As follows from our 
results, the effective interaction in such nonlinear waveg- 
uides is nonlocal, and the nonlinear localised modes are 
described by a nontrivial generalisation of nonlinear lat- 
tice models with long-range coupling and nonlocal non- 
linearity. 



II. MODEL 

We consider a 2D photonic crystal created by a square 
lattice of parallel, infinitely long dielectric columns (or 
rods) in air. The system is characterized by the dielec- 
tric constant e(x) — e(xi, x%), and it is assumed that the 
rods are parallel to the x% axis. The evolution of the TM- 
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polarisecH light [with the electric field having the struc- 
ture E = (0,0, E)], propagating in the {x\, a^-plane, is 
governed by the scalar wave equation 



V 2 E(x,t)-\df[e(x)E]=0, 



(1) 



d 2 



where V 2 

consider the stationary solutions 



d 2 2 . For monochromatic light, we 



E(x,t) = e~ iuJt E{x\u) , 

and the equation of motion ([!]) reduces to the simple 
eigenvalue problem 



E(x\uj) = 



(2) 



This eigenvalue problem can be easily solved (e.g., by the 
plane waves methodliS) in the case of a perfect photonic 
crystal, for which the dielectric constant e(x) = e pc (x) is 
a periodic function 

£pc(*E "h = ^pc(^) : (3) 

where i and j are arbitrary integer, and 

s i:j = ia 1 +ja 2 (4) 

is a linear combination of the primitive lattice vectors a\ 
and a 2 of the 2D photonic crystal. 
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FIG. 1. The band-gap structure of the photonic crystal 
consisting of a square lattice of dielectric rods with ro = 0.18a 
and eo = 11.56 (the band gaps are shaded grey). The top 
right inset shows a cross-sectional view of the 2D photonic 
crystal. The bottom right inset shows the corresponding Bril- 
louin zone, with the irreducible zone shaded grey. 

For dcfinitcness, we consider the 2D photonic crystal 
earlier analysed (in the linear limit) in Refs. [l4],[l5], i.e. we 
assume that the rods are identical and cylindrical, with 
radius ro = 0.18a and dielectric constant eo = 11-56. 
The rods form a perfect square lattice with the distance 
a between two neighbouring rods, i.e. a± = ax\ and 
Gt2 = ax2- The frequency band structure for this type 
of 2D photonic crystal, and for the selected polarisations 



of the electric field, is shown in Fig. [j]. As follows from 
the structure of the frequency spectrum, there exists a 
large (38%) band gap that extends from the lower cut-off 
frequency, u> — 0.302 x 27rc/a, to the upper band-gap fre- 
quency, lu — 0.443 x 2irc/a. Since the characteristics of a 
PBG material remain unchanged under rescaling, we can 
assume that this gap is created either in the infra-red or 
visible regions of the spectrum. For example, if we choose 
the lattice constant to be a — 0.58 /im, the wavelength 
corresponding to the mid-gap frequency will be 1.55 /im. 

The TM-polarised light cannot propagate through the 
photonic crystal if its frequency falls inside the band gap. 
But one can excite guided modes inside the forbidden fre- 
quency gap by introducing some interfaces, waveguides, 
or defects. Here, we consider waveguides created by a 
row of identical defects with a Kerr-type nonlinear re- 
sponse. These defect-induced waveguides possess trans- 
lational symmetry, and the corresponding guided modes 
can be characterized by the reciprocal space wave vector 
k directed along the waveguide. Such a guided mode has 
a periodical profile inside the waveguide, and it decays 
exponentially outside it. 

Linear photonic- crystal waveguides created by remov- 
ing a row of dielectric rods have been recently investi- 
gated numericallyli-3li3 and experimentally.ll3 In particu- 
lar, highly efficient transmission of light, even in the case 
of a bent waveguide, has been demonstrated. 



In the present paper, in contrast to Refs. |14| - |16| where 
only linear waveguides were considered, we study the 
properties of nonlinear waveguides created by inserting 
an additional row of rods fabricated from a Kerr-type 
nonlinear material characterized by a third-order nonlin- 
ear susceptibility with the linear dielectric constant e<j. 
For definiteness, we assume that e^ = eo = 11.56. As 
we show below, changing the radius of these defect 
rods and their location within the crystal, we can create 
waveguides with quite different properties. 



III. EFFECTIVE DISCRETE EQUATIONS 

Writing the dielectric constant e(x) as a sum of peri- 
odic and defect-induced terms, i.e. 

e(x) = e pc (x) + Se{x \ E) , 

we can present Eq. (||) as follows, 

V 2 + (^) 2 e pc (a;) E(x\uj) 

= -(-) 2 5e(x\E)E(x\v) . (5) 
Equation (||) can also be written in the integral form 
E(x | oj) = (^) y d 2 y G(x, y \ u) Se(y \ E) E(y \ u) , 

(6) 
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where G(x, y\u>) is the Green function which is defined, 
in a standard way, as a solution of the equation 



V 2 



.(as) 



G{x,y\uj) = -6(x - y) 



with, accordingly to Eq. (||), periodic coefficients. The 
properties of the Green function and the numerical meth- 
ods for itSpCalculation have been already described in the 
literature Here, we notice that the Green function 
of a perfect 2D photonic crystal is symmetric, i.e. 

G(x,y\u) = G(y,x\u) 

and periodic, i.e. 

G(x + Sij,y + s i:i \u) = G(x,y\w) , 

where is defined by Eq. (||). 

Let us consider a row of nonlinear defect rods embed- 
ded into the crystal along a selected direction. To de- 
scribe such a row, we should define the rods positions 
along Sij with some specific values of i and j. For exam- 
ple, let us first assume that the defect rods are located at 
the points x m = Xq + msij. In this case, the correction 
to the dielectric constant is 



5e(x) = {e d + \E(x\iu)\ 2 }J2Hx-x m ) 



(7) 



where 



6(x) 



for \x\ < rd 
for \x\ > r d 



Assuming that the radius of the rods, rd, is sufficiently 
small (so that the electric field E(x | ui) is almost con- 
stant inside the defect rods), we substitute Eq. (0) into 
Eq. dJ) and, averaging over of the cross-section of the 
rods, derive an approximate discrete nonlinear equation 
for the electric field 



E,, 



E 



Jn-m{u)(ed + \E m \ 2 )E„ 



where 



J„(cj) = (^j J d 2 y G(x , x n + y\ uj) 



(8) 



(9) 



This type of discrete nonlinear equation for photonic 
crystals has been earlier introduced by McGurnliia, for 
the special case of nonlinear impurities embedded in the 
linear rods. However, the analytical approach developed 
by McGurn for that model did not take into account the 
field distribution via the explicit dependence of the cou- 
pling coefficients J n {u>) and, as a result, the equation 
(gj) was not solved exactly. Moreover, the analysis of 
Rcf. [l^ was based on the nearest-neighbour approxima- 
tion where the coupling coefficients are approximated as 
Jn = Jofin,o + JiS n ,±i with constant Jq and J%. 



G(x , x +y I co) 




FIG. 2. The Green function G(xo, xo + y | w) for xo = ai/2 
and u — 0.33 x 2-irc/a. 

In a sharp contrast, in the present paper we provide a 
systematic numerical analysis of different types of non- 
linear localised modes in the framework of a complete 
model. In particular, we reveal that the approximation of 
the nearest-neighbour interaction is very crude in many 
of the cases analysed. Since the effective coupling coeffi- 
cients are defined by the Green function, this can be seen 
directly from Fig. |] that shows a typical spatial profile 
of the Green function which, in general, characterises a 
long-range interaction, very typical for photonic crystal 
waveguides. As a consequence of that .the coupling co- 
efficients \J n (u})\ calculated from Eq. decrease expo- 
nentially with the site number n, and in the asymptotic 
region they can be presented as follows 



Jq(cj) , for n = , 

J*(w)e- a ^)l"l , for \n\ > 1 , 



where the characteristic decay rate a(ui) can be as small 
as 0.85, depending on the values of w, Xq, ay, and rd, 
and it can be even smaller for other types of photonic 
crystals. 

This result allows us to draw an analogy with a class of 
the nonlinear Schrodinger (NLS) equations that describe 
nonlinear excitations in quasi-one-dimensional molecular 
chains with long-range (e.g. dipole-dipole) interactioH. 
between the particles and local on-site nonlinearities.lia 
For such systems, it was shown that the effect of nonlo- 
cal interparticle interaction introduces some new features 
in the properties of existence and stability of nonlinear 
localised modes. Moreover, for our model the coupling 
coefficients J n (w) can be either non-staggered and mono- 
tonically decaying, i.e. J n {oj) = |J„(w)|, or staggered and 
oscillating from site to site, i.e. J n (uj) = (— l) n \ J n {u>)\. 
We can therefore expect that effective nonlocality in both 
linear and nonlinear terms of Eq. (||) will bring a num- 
ber of new features in the properties of nonlinear localised 
modes. 
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IV. EXAMPLES OF NONLINEAR MODES 



As can be seen from the structure of the example Green 
function, presented in Fig. the case of monotonically 
varying J n {w) can be obtained by locating the defect rods 
at the points Xo = ai/2, along the straight line in the Sqi 
direction. In this case, the frequency of a linear guided 
mode, that can be excited in such a waveguide, takes the 
minimum value at k = (see Fig. |J), and the correspond- 
ing nonlinear mode is expected to be non-staggered. 
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FIG. 3. Dispersion relation for the photonic crystal waveg- 
uide shown in the inset (eo = ed = 11.56, ro = 0.18a, 
rd = 0.10a). The grey areas are the projected band struc- 
ture of the perfect 2D photonic crystal. The frequencies 
at the indicated points are: lua = 0.378 x 2irc/a and 
ujb = 0.412 x 27rc/a. 
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FIG. 4. Intensity Q(lu) of the nonlinear mode excited in the 
photonic crystal waveguide shown in Fig. |j| The right inset 
gives the dependence J n (ui) calculated at u) = 0.37 x 2-nc/a. 
The left inset presents the profile of the corresponding non- 
linear localised mode. 



the guiding properties of a channel created by defect rods, 
and in the direction of the waveguide, due to the self- 
trapping effect. Such nonlinear modes exist with fre- 
quencies below the frequency of the linear guided mode 
of the waveguide, i.e. below the frequency lua in Fig- 
and are indeed non-staggered, with the bell-shaped pro- 
file along the waveguide direction shown in the left inset 
of Fig. a 
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FIG. 5. Dispersion relation for the photonic crystal waveg- 
uide shown in the inset (to = £d = 11.56, ro = 0.18a, 
r<j — 0.10a). The grey areas are the projected band struc- 
ture of the perfect 2D photonic crystal. The frequencies 
at the indicated points are: lua = 0.346 x 2irc/a and 
uj B = 0.440 x 27rc/a. 




0.32 0.33 0.34 0.35 

Frequency (ooa/27tc) 

FIG. 6. Intensity Q(lu) of the nonlinear mode excited in 
the photonic crystal waveguide shown in Fig. |f| The solid 
curve corresponds to the case of nonlinear rods in a linear 
photonic crystal, whereas the dashed curve is the same de- 
pendence for the case of a nonlinear photonic crystal. The 
right inset shows the behaviour of the coupling coefficients 
J„(lu) for n > 1 (J = 0.045) at u = 0.33 x 27rc/a. The left 
inset shows the profile of the corresponding nonlinear mode. 



We have solved Eq. (||) numerically and found that 
nonlincarity can lead to the existence of a new type of 
guided modes which are localised in both directions, i.e. 
in the direction perpendicular to the waveguide, due to 



The 2D nonlinear modes localised in both dimensions 
can be characterized by the mode intensity which we de- 
fine, by analogy with the NLS equation, as 
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Q = Yl \ E " 



(10) 



This intensity is closely related to the energy of the elec- 
tric field in the 2D photonic crystal accumulated in the 
nonlinear mode. In Fig. |J we plot the dependence of Q 
on frequency, for the waveguide geometry shown in Fig. 
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FIG. 7. Dispersion relation for the photonic crystal waveg- 
uide shown in the inset (eo = = 11.56, ro = 0.18a, 
rd = 0.10a). The grey areas are the projected band struc- 
ture of the perfect crystal. The frequencies at the points in- 
dicated are: lua = 0.352 x 2irc/a, lob = 0.371 x 2nc/a, and 
ojc = 0.376 x 2-Kc/a (at k = 0.217 x 2-ir/a). 
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FIG. 8. Intensity Q(lo) of the nonlinear mode excited in 
the photonic crystal waveguide shown in Fig. (?]. The right 
inset shows the behaviour of the coupling coefficients J n (co) 
for n > 1 (Jo = 0.068) at cj = 0.345 x 27rc/a. The left inset 
shows the profile of the corresponding nonlinear mode. 

As can be seen from the example of the Green function 
shown in Fig. the case of staggered coupling coeffi- 
cients J n {<J) can be obtained by locating the defect rods 
at the points x$ — ai/2, along the straight line in the 
Sio direction. In this case, the frequency dependence of 
the linear guided mode of the waveguide takes the mini- 
mum at k = ir/a (see Fig. ||). Accordingly, the nonlinear 



guided mode localised along the direction of the waveg- 
uide is expected to exist with the frequency below the 
lowest frequency to a of the linear guided mode, with a 
staggered profile. The longitudinal profile of such a 2D 
nonlinear localised mode is shown in the left inset in Fig. 
^, together with the dependence of the mode intensity Q 
on the frequency (solid curve), which in this case is again 
monotonia 
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FIG. 9. Contour plot of the the nonlinear localised mode 
that corresponds to the longitudinal cross-section shown in 
the left inset of Fig. |. 

The results presented above are obtained for linear 
photonic crystals with nonlinear waveguides created by 
a row of defect rods. However, we have carried out the 
same analysis for the general case of a nonlinear photonic 
crystal that is created by rods of different size but made 
of the same nonlinear material. Importantly, we have 
found very small difference in all the results for relatively 
weak nonlinearities. In particular, for the photonic crys- 
tal waveguide shown in Fig. ^, the results for linear and 
nonlinear photonic crystals are very close. Indeed, for 
the mode intensity Q the results corresponding to a non- 
linear photonic crystal are shown in Fig. ^ by a dashed 
curve, and for Q < 20 this curve almost coincides with 
the solid curve corresponding to the case of a nonlinear 
waveguide embedded into a linear photonic crystal. 

Let us now consider the waveguide created by a row of 
defects which are located at the points Xq = (ax +a,2)/2, 
along a straight line in either the Sio or Sqi directions. 
The results for this case are presented in Figs. f^]. The 
coupling coefficients | J n \ are described by a slowly decay- 
ing function of the site number n, so that the effective 
interaction decays on scales much larger than those in the 
cases considered previously. Similar to-jLhe NLS models 
with long-range dispersive interactionsEfrtj, for this type 
of nonlinear photonic crystal waveguide we find a non- 
monotonic behaviour of the mode intensity Q(u>) and, as 
a result, multi-valued dependence of the invariant Q(lo) 
for uj < 0.347 x 2nc/a. Similar to the results earlier 
obtained for the nonlocal NLS modelslia, we can expect 
here that nonlinear localised modes corresponding, in our 
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notations, to the positive slope of the derivative dQ/dui 
are unstable and will eventually decay or transform into 
modes of higher or lower frequency. Such a phenomenon 
is known as bistability, and in this problem it occurs as 
a direct manifestation of the nonlocality of the effective 
(linear and nonlinear) interaction between the defect rod 
sites. 



V. CONCLUSIONS 

Exploration of nonlinear properties of PBG materials 
is a new direction of research, and it may open up a new 
class of applications of photonic crystals for all-optical 
signal processing and switching, allowing an effective way 
to create tunable band-gap structures operating entirely 
with light. Nonlinear photonic crystals, and nonlinear 
waveguides embedded into photonic structures with pe- 
riodically modulated dielectric constant, create an ideal 
environment for the generation and observation of non- 
linear localised modes. 

In the present paper, we have developed a consistent 
theory of nonlinear localised modes which can be excited 
in photonic crystal waveguides of different geometry. For 
several geometries of 2D waveguides, we have demon- 
strated that such modes are described by a new type of 
nonlinear lattice models that include long-range interac- 
tion and effectively nonlocal nonlinear response. It is 
expected that the general features of nonlinear guided 
modes described here will be preserved in other types of 
photonic crystal waveguides. Our approach and results 
can also be useful to develop the theory of nonlinear two- 
frequency parametric localised modes in the recently fab- 
ricated 2D photonic-crystals with the second-order non- 
linear susceptibilityEj. Additionally, similar types of non- 
linear—localised modes are expected in photonic crystal 
fibersEll consisting of a periodic air-hole lattice that runs 
along the length of the fiber, provided the fiber core is 
made of a highly nonlinear material (see, e.g., Ref. |22|). 
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